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INTRODUCTION 

For many problems it is reasonable to expect that the minimum time solution 
is also the minimum fuel solution. However, if one allows the propulsion system 
to be turned off and back on, it is clear that these two solutions may differ. 
In general, high thrust transfers resemble the well-known impulsive transfers 
where the burn arcs are of very short duration. The low and medium thrust 
transfers differ in that their thrust acceleration levels (10 4 earth g’s for low 
thrust and 10 — 2 earth g’s for medium thrust accelerations) yield longer burn 
arcs which will require more revolutions, thus making the low thrust transfer 
computational intensive. Here, we consider optimal low and medium thrust 
orbital transfers. 


APPROACH 

The formulation of the problem follows that of Horsewood et. al. (1), where 
the state of the spacecraft is given in terms of the slowly varying equinoctial 
elements, which are expressed in terms of the classical elements a,e,i,w, Q as 

T 

z = (a, e sin (u> + fl), e cos (u> + fl), tan (i/2) sin fl, tan (i/2) cos f l) 

The spacecraft mass, m, is also a state variable. In addition, the position of the 

spacecraft within an orbit is given by the eccentric longitude, F = E + u + fl. 

The equations of motion are 

• 2P . • 2P m 

z = — Mu, m = — — [ij 

me c~ 

where M = || isa5x3 matrix calculated by treating F as an independent 
parameter, such that the variation of F with respect to the other state variables 
is zero; P is the power due to the thrusters, c is the exhaust velocity and u is 
the unit vector in the direction of thrust. 


It follows from the well-known maximum principle that the fuel optimal tra- 
jectory from a given state to some desired final state is found by thrusting at 
every point along the trajectory in the direction which maximizes the Hamil- 
tonian function H. H can be written in terms of the state variables and their 
corresponding costate (adjoint) variables A z and A m as 


H = A z T z -I- A m rn 

me \ c ) 


[ 2 ] 
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where the costate variables satisfy a first order linear system of ordinary dif- 
ferential equations. Now, to maximize H we need only thrust in the direction 
given by u = Note that if the quantity in parenthesis in [2], call it <7, is 

negative, then H is negative and hence H is maximized by letting P = 0, which 
amounts to turning the propulsion system off , i.e. “coasting”; on the other 
hand, if a is positive, then H is maximized by letting P take on its maximum 
value, i.e. “thrusting”. 

Because of the many orbit revolutions of a long duration transfer, the 
intensive computations can be reduced by an averaging technique. We can 
compute an “averaged” Hamiltonian function by holding the state and costate 
variables constant over an orbital period of duration r, i.e. we assume Keplerian 
motion, and integrating the actual Hamiltonian function as follows: 

T 

H = - J 

0 

2 7T 

= J H (z,\ z ,m,\ m , F)s(z, F)dF 

0 

where s(z, F) = ^ . The "averaged” equations of motion can now be computed 

using this “averaged Hamiltonian, and since the Hamiltonian and its derivatives 
are zero during coast phases, we need only integrate these equations over the 
predetermined thrust intervals. 

NUMERICAL APPROACH 

A FORTRAN program MINFUEL has been written to compute this 
minimum-fuel solution. The basic algorithm is: 

(1) start iteration of the unknown initial costate values by making a first 

guess for the simplified problem of a circle- to-circle coplanar transfer, an initial 
costate guess can be found by using the fact that it is fuel-optimal to thrust in 
the direction of motion, or tangential direction. This known direction can be 
used to help guess the initial costate values since u = Given any state, 

this optimal thrust direction can be calculated in terms of the costate values by 
running the program XFORM to compute the transpose of the matrix M and 
the velocity vector. 

(2) compute a trajectory corresponding to an initial costate guess — one calls 
the Runge-Kutta routine to integrate the averaged equations of motion; the 
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switch function can be expressed in terms of the fast variable F around an orbit 

and the zeroes of the switch function a = |a^MM t A z 2 are found, which 

defines the thrusting subintervals of [0, 2x] for integration. Again for the circle- 
to-circle transfer we know that initially the switch function should be positive. 

(3) compare the computed final state with the specified final state and the 
proper transversality conditions, such as \ mj = 1. The Newton method used for 
this iteration is very sensitive to the initial costate guess. 

The medium thrust problem was investigated through the use of the Flight 
Mechanics Branch’s program SCOOT (Simplex Computation of Optimal Orbital 
Transfers). SCOOT was known to converge to a quasi-optimum solution for high 
thrust levels and short bum times. Although low thrust transfers require long 
burn times, SCOOT was thought to be a viable option for guessing optimal low 
thrust trajectories by using the final output from the converged higher thrust 
cases as input to MINFUEL. After many runs of SCOOT, the fact that it 
calculates only a local optimal and is thus sensitive to coast and burn times 
guessed, makes it very difficult for comparisons of medium and low thrust 
transfers. 

CONCLUSIONS 

The results of SCOOT have indicated that lowering the thrust level “can” 
increase the amount of final propellant; increasing the ISP can increase the final 
propellant as well. The output of SCOOT, in particular the optimal thrust 
directions, from medium thrust level solutions can be used to help calculate an 
initial guess for low thrust transfers. The program MINFUEL, once properly 
debugged, should prove very capable of handling minimum-fuel low thrust 
transfers. 

REFERENCES 

(1) Horsewood, J.L., Suskin, M.A., Pines, S., Moon Trajectory Computa- 
tional Capability Development. Technical Report 90-51, NASA Lewis Research 
Center, July 1990. 

(2) Jasper, T.P., Low-Thrust Trajectory Analysis for the Geosynchronous 
Mission. AIAA 10th Electric Propulsion Conference, 1973. 

(3) Sackett, L.L., Malchow, H.L., Edelbaum, T.N., Solar Electric Geocentric 
Transfer with Attitude Constraints: Analysis. Technical Report NASA CR- 
134927, The Charles Stark Draper Laboratory, Inc. August 1975. 


XI-3 



